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We investigate the magnetic properties of quasi one-dimensional quantum spin-S antiferromag- 
nets. We use a combination of analytical and numerical techniques to study the presence of plateaux 
in the magnetization curve. The analytical technique consists in a path integral formulation in terms 
of coherent states. This technique can be extended to the presence of doping and has the advantage 
of a much better control for large spins than the usual bosonization technique. We discuss the ap- 
pearance of doping-dependent plateaux in the magnetization curves for spin-5 chains and ladders. 
The analytical results are complemented by a Density matrix Renormalization Group (DMRG) 
study for a trimerized spin- 1/2 and anisotropic spin-3/2 doped chains. 
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I. INTRODUCTION 

One-dimensional and quasi one-dimensional quantum 
spin systems have been a subject of close attention 
in condensed matter physics. Over the past decades, 
much progress has emerged from the studies of such low- 
dimensional systems, the Haldane gapi and the magneti- 
zation plateauJs^Ti^ being some of the best-known exam- 
ples. In particular, spin chains have allowed to study 
quantum magnetism in the simplest tractable setting, 
while investigations of spin ladders have permitted first 
steps in the study of crossover from one to two dimen- 
sions. 

The present work is devoted to magnetization plateaux 
in quantum antiferromagnets: the magnetization remain- 
ing spectacularly constant in a finite interval of external 
magnetic field. This phenomenon has been found in a 
great variety of systems from spin ladders^ and p-merizcd 
chains^ to frustrated higher dimensional systems i^. 

For a quantum spin-S* chain, the necessary condition 
for the occurrence of a magnetization plateau has been 
established by Oshikawa, Yamanaka and AffleckSi: 



N{S-m) e Z, 



(1) 



where m is the magnetization per site, N is the number 
of spins per unit cell, and Z is the set of all integers. 

This condition restricts the plateau magnetization m 
to rational values. However, it has been argued that 
plateaux may also appear at an irrational m, as a re- 
sult of either quenched disorder— or doping with itinerant 
carriers^"— i^^d^ . In the latter case, doping may allow to 
reduce the m in a controlled way, thus making the plateau 
more easily accessible to experiments in lower magnetic 
fields. 

Unfortunately, the analytical methods, used so far to 
study magnetization plateaux in spin chains and lad- 
ders, have been limited to bosonization (and thus, effec- 
tively, to one-dimensional spin-1/2 systems) - and to the 
bond operator technique, intrinsically restricted to spin- 
1/2 systems^. However, recently, Tanaka, Totsuka and 



Hui^ (TTH) have arrived at the necessary condition of 
Ref.i using the Haldane's path integral approach. This 
approach is applicable regardless of the dimensionality 
or of the value S of the spin - and, below, we explore 
how much progress it may afford us in understanding the 
magnetization plateaux in various systems from spin-1/2 
chains to higer-spin systems. This approach open new 
perspectives in understanding the physics of magnetiza- 
tion plateaux in higher spin and higher dimensional sys- 
tems such as rare earth tetraboridea^^. 

Our goal in this paper is to test and extend the ap- 
proach developed by Tanaka et alt^. First, we test the 
technique on several concrete small-spin N-leg ladders 
and p-merized chains, previously studied with the help 
of bosonization^i^. Then we extend the technique by 
combining it with ideas due to R. Shankar— , to account 
for hole doping in spin chains at non-zero magnetiza- 
tion. Unfortunately, accounting for hole doping in zero 
field has proven to be problematic for technical reasons 
which we arc not going to discuss here. However, at 
a non-zero average magnetization per site the technical 
difficulties are lifted, which allows us to generalize the 
TTH approach to the case of doping. After this, we use 
doped Hubbard and t — J chains and ladders as a test- 
ing ground, and re-derive some of the results previously 
found for these systems^"— lii. Then, we turn to doped 
higher-spin systems, which tend to pose a problem for 
the bosonization approach. 

In Section [U we illustrate the key points of Ref. [l3 bv 
studying magnetization plateaux in an anisotropic spin 
chain with the help of coherent-state path integral tech- 
nique. In Section Um we extend this to a dimerized chain, 
and to a two-leg spin ladder - and then, in section IIV| 
we generalize the above to n-leg ladders and p-merized 
chains, and find that they satisfy the necessary condition 
of the Eq. (P). 

In Section |Vl we generalize the plateau condition in 
the Eq. ([1]) to account for the doping dependence. In a 
spin-5 system at a small hole density 5^1, we find the 



plateau condition to read 
5 
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The Section IVII presents an application of the formal- 
ism to a doped spin-3/2 chain, and compares the results 
with those obtained numerically by DMRG. 

In Section IVIH we discuss the plateau condition for 
a doped n-leg ladders and p-merized chains, with and 
without doping. 

In Section I Villi we present the results for a trimerized 
chain, and compare with DMRG results. 

Finally, in Section HXl we present the conclusions, pos- 
sible implications and extensions of the present approach 
to higher dimensions. 



II. ANISOTROPIC SPIN CHAIN IN A 
MAGNETIC FIELD 

In this Section, we study the spin-S* nearest-neighbor 
antiferromagnetic (AF) chain with easy-plane single-ion 
anisotropy, subject to a transverse magnetic field: 

i j 3 

where J is positive, and the magnetic field h points in 
the z direction, as shown in the Fig. [TJ 

Following the Ref. [13, we analyze the system using the 
coherent-state path integral description due to Haldanoi. 
The resulting effective action comprises two terms: the 
first one is the coherent-state expectation value of the 
Hamiltonian, and the second one, dubbed the Berry 
phase term, corresponds to the solid angle swept by the 
spins in their imaginary time evolution. 

In order to obtain an effective theory, first, we iden- 
tify the classical ground-state configuration and the low- 
energy modes above it. Partially polarized by the mag- 
netic field, the spins form a canted texture, that we 
parametrize as Sj = Sfij^ with fij being a unit vector 



with staggered XY components 



Sj = {Ssindj cos( 



)j , S sin 0j 



smc 



<t)j,Scosej) . (4) 



We parametrize the fluctuations around the above 
canted state as per 



TT 

—Xi 

a 



(5) 



where is the classical ground state solution cos^o 
h J ^ _ with a the lattice constant. 



25(2,/- 



Expanding up to the second order in the 56, we can 



write S, 



Sj zt iSj as a function of > 



„^ ^ (xj) and (l){xj) 

and, using these fluctuation fields, we write an effective 
theory. If we calculate the Poisson Brackets {S^, S^}^,s0, 
we obtain in{S'^S'±}0,5e = -5(sin6lo-(56'cos6'o) (±^5*^) 
then is straightforward to see that defining 



all{xj) = —5 



^^(xj) sin^o + cos 00 



(6) 
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FIG. 1. Classical configuration for a spin chain in the presence 
of a magnetic field. 



as the conjugate field of we have the correct commuta- 
tion relations for the spin operators 



S^ ~ Scos0o + an{xj) 



a2 52 



S sin 6*0 
I 



5* sin do 



(7) 

n(x-j) (8) 
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where m = S'cos^o- 

Following the Ref. [ij, one arrives at the low-energy 
continuous effective action, corresponding to the Hamil- 
tonian ([3]): 



S^J dxdr {a{2J + D)I1^ -i{dr<l))n 



-a{S' ™2^,« ^^2 
2 ^ 



(9) 



The last term in right hand side of the Ec^. ^ arises from 
the Berry phase of the individual spins in the Eq. 
After gaussian integration over the field 11, the action 
(l9l) takes the form 



S-- 



dxdr 



2 



with 



1 



2a (2 J + D) 



^ Ja{S'^ 



(10) 



The last term in (|10p counts the winding number of 
the space-time history of the field 0, defined on a cover- 
ing space of a circle. In order to understand the con- 
sequences of the topological term, it is convenient to 
apply a standard duality transformation to the action. 
First, the phase field cj) is written as 4> = 4>y + (pt, where 
(j)v is a fixed field configuration containing all the vor- 
tices {df^di, — d,ydf_i)(j)v ^ 0, and (j)t contains the fluctu- 
ating vortex-free part. Next, we introduce the Hubbard- 
Stratonovich auxiliary vector field = {Jt^Jx), and 
integrating by parts we obtain 



S 



dxdr 



i\ Jr 



S — m 



1 J2 1 J2 
+ ^J.(a.0.) + -^+-^ 



(12) 



dr Jt 



s- 



dxJx 
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Defining J^, = + 5, 



-) the vorticity-free part can 



be eliminated with the constraint 9^ = 0. Then, the 
action reads 



5" -771 



The constraint 5^ = is solved in one dimension (ID) 
in terms of an auxiliary field by 



(13) 



where x is vorticity-free. Integrating the first term in 
the action by parts, and using the substitution x ~^ X + 
{^~a^) '''' obtain 



+ iB [x + I — — — I X 



(14) 



where B ~ t^vd^dycj)^ — 2'n:qj5{T—fj)5{x—Xj). Here, 
Tj and Xj are the time and space coordinates of the j-th 
vortex event and qj is the vorticity. Upon summation 
over the vortex configurations in the partition function, 
and then rescaling the time variable, the action is brought 
to the form 



S — dxdr 



+ Ai cos 27r X + 



5-777 



(15) 



This effective action is identical to the one obtained in 
the Ref. [2 via bosonization. The cosine term is com- 
mensurate for S — m € Z. The magnetic excitations arc 
gapless for all the values of m that satisfy S — rn ^ Z. 
Then, when the condition 



5 — TTi € Z 



(16) 



is satisfied, a plateau can occur for some range of pa- 
rameters if the scaling dimension of the perturbation is 
small enough. More precisely, the presence of the cosine 
operator in the Eq. is not enough to produce a gap 
in the spectrum. The stiffness of the x held must be such 
that this operator has scaling d dimension smaller than 
2. Since d is a decreasing function of D, we expect a 
plateau for large enough values of this parameter. 

The condition is the well-known result of Os- 

hikawa, Yamanaka and Affleck^, which generalizes the 
argument due to Lieb, Schulz and Mattisii. 

A more general situation arises when S — m = p/r, 
with p,r € Z. In this case, the vortices whose winding 
numbers are integer multiples of r are free of destructive 
interference. In this case, only vortices with vorticity ±r 
are able to condense. In the effective action, this is sig- 
nalled by the rth harmonic of the cosine term in the har- 
monic expansion being commensurate. If relevant, this 
cosine operator gives rise to a gapped r-fold degenerate 
ground state and fractionalized excitations^^. The best- 
known example of such a scenario is the AF chain with 



a strong enough second-neighbor AF couplingi^. The 
system is then in a gapped phase with two degenerate 
spontaneously dimerized ground states; fractional exci- 
tations are spinous - domain walls between these two 
ground states. In what follows, we concentrate on the 
case S — m e Z. 

At this point, we would like to make an observation, 
which will prove useful in the Section|Vl when comparing 
the plateau condition for doped systems with that at zero 
doping: According to the Eq. pl}|) . for each contribution 
to the partition function, the Berry phase term provides 
a weight factor 

giiS-m}J2j J dTdr4>j ^ 

where the integral J dTdrfj^j = 2Tmj is an integer times 
2tt, and the subscript j labels the spins in the chain. 
Inverting the sign in front of the S changes the imagi- 
nary exponent by AnrijS, which is an integer times 2tt 
both for integer and half-integer S. Therefore, the nec- 
essary condition (|16|1 for the magnetization plateau may 
be equivalently presented as 



(5±m) e Z. 



(17) 



In the following, we would like to extend the approach 
above to other one-dimensional spin systems. In the next 
section, we discuss the cases of a two-leg ladder and a 
dimerized spin chain. Then, in the Section [Vj we study 
the magnetization plateaux in the presence of doping. 



III. TWO-LEG LADDERS AND DIMERIZED 
CHAINS 

A. Two-leg spin ladder 

The formalism above can be used to study more in- 
volved spin models, such as spin ladders, that interpolate 
between one and two dimensions. For spin 1/2 experi- 
mental evidence of zero magnetization plateaux has been 
reported for instance in (C5Hi2N)2CuBr4i^. In this sec- 
tion we study two-leg ladders with single-ion anisotropy 
in a magnetic field, and in the next section we discuss 
extensions to iV-leg ladders and p-merized chains. 

Consider the following Hamiltonian 



H : 



Si 



J+1 



+ S2,, ■ S: 



2,3 + 1) 



(18) 

- h{{Sl,) + {SI,))} , 

where J,, is the antiferromagnetic coupling along the 
chain, and J± is the inter-chain coupling. The onsite 
anisotropy term is added only for completeness and we 
can recover the isotropic case by taking the D ^ limit. 

In the S ^ oo limit we can write the energy of the 
system as a function of the angle 9 (see Fig. [2]), with the 
minimum at 

^°^'°^ 25(2J„+V+i^) - ^''^ 



As before, we parametrize the fluctuations in terms of 
the fields Sda and 0q , where a = 1,2 is the chain in- 
dex. Equations ([7]) and ([5]) have the same form in 
each chain. Using this parametrization in the Hamilto- 
nian, retaining terms up to second order in the fields and 
switching to the path integral language, we obtain the 
action 

S = J rfsdr|ya(S'-m')((9,0i)' + (9,(?!.2)') 



~^ 2 \ a 



+ J dxdr |a y2J„ + J 



+ D] (n? + n^) 



52 - m2 



(20) 



Integrating over the H-ficlds and using the transforma- 
tion Lp = Ucf) with 



if 



U 
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(21) 



we can write S = Ss + Sa, with 

5. = 



drdx i iif(«)(a,^,)2 + (5,(^3 )2 



i 2 
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and 



"5'a 



dxdr Uxi-Hd^^a)' + ^Ki^Hdr^aY 



2J_L 
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where 



4J„ + Jx ( 1 + 3 



4J„ + Jx 3 + 
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The action corresponding to the antisymmetric field 
ipa contains a mass term. In order to obtain an effective 
theory the field (fa can be evaluated in the saddle point 
solution ipa = 0. 

In the symmetric action Ss, we use the duality trans- 
formation of the section [U presenting the field ips as 
a sum of the vortex component (ps.v and the vortex- 
free component ips,t- Upon rescaling the time as per 

T \J K'f^/Kx^'^ T, and following the standard steps 
we obtain 
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FIG. 2. Classical configuration for a spin ladder in the pres- 
ence of a magnetic field. 



Now, the cosine term is commensurate if 
2(5" ±771) e z 



(23) 



Note that, even though our approach does not apply 
directly to the ttt, = case, this conditions is nevertheless 
consistent with the well known m = plateau for S = 
1/2. 

Another important comment is due here. The factor 2 
in front of the Berry term is not an artifact of the trans- 
formation (PT|) because of the periodicity of the fields. 
Each one of the fields 0i and 4>2 satisfy 



(j)i{x + L) = 01 (.t) + 27r7ii 

(/)2(.T + L) = (f)2{x) + 27r7i2 



where ni and ni are integers. Then the antisymmetric 
combination satisfy (j)i — 4>2 = 27r(77i — 712). The saddle 
point solution (pa = Q implies ni = n2 and then the 
sum (f>i + 02 has periodicity Anni. The factor i in the 
definition of ips is necessary for the correct periodicity. 



B. Dimerized spin chain 

The dimerized spin chain can be studied in a similar 
way. We begin with the Hamiltonian 



H = ^ {(J + 5J)Si.j ■ S2,, + (J - 5J)S2,j ■ Si,,+i 
+ D [{Sl^f + [Sl^f] - h [St, + S^,] } , (24) 



S — dxdr ■ 



2\/Ki'^Ki'^ 
+ Ai cos ( 27r ( X + 2 ( ^—H!: I X 



(22) 



and use the Eqs. ([6]), ([7]) and ([8]) for each sublattice. 
We then switch to the path integral formalism, take the 
continuum limit and drop the constant terms to obtain 
the following action: 



5 



S ~ clxdr < J 



+ -(J ~ 6J)iS' - m') + (9,^2)'] 

+ i( J - sj){s^ - m^){^^ - 02)(a,(/.i + d,<P2) 



a J 



g2 _ jyi2 



2aJ 1 



D]{ni + ni) 



U1U2 



(25) 



Integrating the H-ficlds and making the substitution (f 
U(j) with U given by pip wc have 



dxdr < 4 J 



^2- 



+ -(J - 5J)(52 - m^) [{d,^,f + id.ipaf] 

' S — 



i2 



(26) 



The field is massive and we can use the saddle point 
solution for it ipa = 0. The action for the symmetric field 
results 



Ss = / drdx { ^K^id^ips)^ + ^kr{dr^s? 
2l [ 1 (drips) 



(27) 



with = a{J ~5J){S^ ~rr?) and = l/[a(2J + D)]. 

The action has the same form as for the two-leg lad- 
der and, repeating the steps described in the preceding 
subsection, we obtain the necessary condition for the for- 
mation of magnetization plateaux 



2{S±m) e Z. 



(28) 



In both cases studied in this Section, the unit cell con- 
tains two spins, and the effective action was written in 
terms of two fields: and ipa- Notice that, of these two, 
only the massless one (i^^) defines the necessary condition 
for the formation of plateaux. 



IV. EXTENSIONS TO N-UEG LADDERS AND 
P-MERIZED CHAINS. 

The arguments, presented in the preceding sections, 
can be easily extended to more complex models. As we 



have shown in the previous section, the magnetization 
processes of a two-leg ladder and of a dimerized chain are 
described by a single effective action of the massless field, 
that has the same form for these two different models. 
Under certain restrictions, this similarity holds for iV-leg 
ladders and A^-merized chains, as well. In the last section 
we have used two fields to describe the low-energy theory 
and finally only one field remains massless in our effective 
action. The extension to the A^-leg ladder is natural and 
follows the same steps as in the two legs case but now 
working with N different fields. The first N —1 fields are 
massive as in the ease of for the two legs ladder. The 
effective action can be written also in terms of the last 
field, which is given hy = ^^'^ ^^'•^ action is 

a straightforward generalization of the of equation 

27r2 



S= / dxdT{ — {{d^^f + {d^^f 

+ Ai COS ( 27r ( X + A*' ( ^ \ x 



(29) 



where /3 depends of the microscopic parameters and the 
commensurability condition is given by 



N{S±m) e Z. 



(30) 



Again, this result contains also what is known for the 
zero magnetic field. The A^-leg ladder at zero field was 
studied by G. Sierrap2£ using the original Haldane's path 
integral approach, and the presence of a spin gap was 
recovered from the condition (l30l) with m ~ 0. 



V. ADDING HOLES TO A SINGLE CHAIN 

There are several examples of hole-doped magnets. For 
instance, doping the S = \ metal oxide Y2BaNi03^ 
with Ca introduces hole carriers in the chains. Other 
examples are given by manganese oxides such as 
Lai-jjCa^jMnOj^. 

Theoretical studies of these compounds generally de- 
part from the double exchange model^^. Using strong 
Hund's rule coupling between the itinerant and localized 
spins, an effective Hamiltonian on a restricted Hilbert 
space is introduced, where creation of a hole on a given 
site replaces the spin S on this site by S ~ 1/2. Such 
a calculation of the effective Hamiltonian generalizes the 
derivation of the t — J model from the Hubbard model for 
S = l/32i. In this section, we study the effect of doping 
on the low energy physics of a Heisenberg chain and the 
corresponding effect on the magnetization curve. 

We begin with a spin chain with one spin per site in 
the presence of a magnetic field. Using the path inte- 
gral formulation, we have shown that the effective action 
is given by the Eq. ([9]) of the Section [III Creation of 
a hole at a given site j corresponds to extracting the 
spin from this site. To account for this, we will simply 
remove the contribution that the jih. spin would have 
made to the action. Let us introduce at each site a hole 
creation operator i/'j , satisfying fermion commutation re- 
lations {V'li V'j} = Without any holes, we have just 



6 



the low-energy action ^ of the Section |TT] for the pure 
spin system. 

Important contributions to the action upon creation 
of a hole are obtained by the change in the interaction 
between the spins, and simply replacing the Berry term 
for the spin-iS* by that for S — 1/2 for the spin with the 
hole and then taking into account the contribution to the 
action due to the hole hopping. We start by considering 
this very hopping term. 

Let rti and ^2 be the classical spins at sites 1 and 
2. According to the arguments of Shankar— , the correct 
hole hopping term from site 1 to site 2 would not be 
simply —tip24'i: such a term would move the hole, but 
not the spin fl. At the microscopic level, the electron is 
transferred by the operator d = X^o- ^1(t^2ct, that moves 
the charge and the spin coherently from 2 to 1. The 
action of this operator on a state with the hole on site 1 



d\0,n2) = 1^2,0) 



(31) 



where o represents a hole. In this language, the correct 
matrix element for the process is —t{D,i\D,2) . Thus, the 
correct hopping term is — t(f2i|f22) V'lV'i: where 



; — cos — + e ^^'^ sm — sm — 



2S 

(32) 



Notice that, in the h — Q case, the classical configuration 
of the spins is antiparalel and the overlap of coherent 
states vanishes. Then, in the zero magnetic field case, the 
hoping amplitude between nearest neighbors is a fiuctu- 
ating variable with zero average, a very difficult problem 
to study. This is the main reason for the original problem 
on hole dopingi^ to concentrate in a model where doping 
between second neighbors was the dominant effect. In 
our case, the non zero magnetization implies a non zero 
overlap between neighboring coherent states allowing for 
a consistent study of first neighbors hoping problem. 

Finally, the hopping of holes is described by the fol- 
I 



lowing Hamiltonian 

Hhopp = -i ^(%|r2j+i) + h.i 



(33) 



If we expand 9 around the classical energy minimum 
value 9(x) = Oq + S9{x), we must retain the terms up 
to order a, since the continuum limit ipj — > y/ail){x) for 
the fermion operators involves an extra factor of a. Using 
the Eq. ([6]) and retaining terms up to the first order in 
a, we find 

^ (^)'''e(2''^n(^)-*'^^(^-")^^*(^))(34) 

then 

Hhop^ty^ (l + 2a—Yl{x) -ia — {S -m)dx(l>{x 



+ h.c. (35) 
Since the sought long-distance physics involves only 
the states near the Fermi surface, we linearize the theory 
around ±kp 



R,3 



-ikp-aj 



(36) 



to obtain Hhopp ~ Hfree + SH, where 

Hfree = ^ -t COS kpa V"]?,-,-)-! V"/? J + i'R^li^RJ+l 



it sin kpa 



i/'flj+iV-flj - 1pR,l^R,3+l (37) 



and 



-2a^{S - ?n) sin {kFa)dx(t>{x) (^V'aj^b.j - '^l^'^ij) | (38) 



Taking the continuum limit -^■ ^/aip{x) and setting 
the Fermi energy to zero, we obtain 



^hopp / dx 



(Til \ -^^ 
— j sin (/cFffl) 



^r{=^) l-id. + ^^^-^d.A ^r{x) - ^l{x) f-id^ + ^^^—^d^4> ) ^l{x) 



AatSm 



cos {kFa)Il{x) (rpl^{x)i>R(x) + i!l^(x)ipL{x)^ j 



(39) 



Since the linearized theory has an infinite number of par- 
ticles in the ground state, we shall introduce normal or- 
dering, to correctly define the theory: 

^t(x)^(x) - 5+ : ^b^j^ix)^PRix) + 4,l{x)Mx) ■■ (40) 



gin with the action without holes, and remove the Berry 
phase term at each hole site. For larger values of S, 
we have more than one electron on each site and a hole 
corresponds to a spin [S — 1/2] impurity in the spin- 5* 
hos t^'^i^^'^^ . In other words, we apply the projection op- 



Now, we account for hole doping. For 5* — 1/2, wc be- 
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erator 



2S' 



(41) 



to the Berry phase terra in the action if the hole was 
created on site j. We must do the same for the spin-spin 
interaction part using the projector. 



2S 



)(1 



25 



(42) 



Note that this is somehow a different scenario than the 
one proposed by Shankar— where an entire spin S would 
hope from one site to another. Our approach is more 
appropriate to cope with the experimental situation we 



mention below. The modification of our approach to han- 
dle Shankar scenario, or any intermediate value of the 
spin hoping between S and 1/2 is straightforward. 

The following step is to linearize the theory near ±kp, 
as we did in the hopping term, to obtain 



1 



2S 
_S_ _ 

25 " 



(43) 



1 

25 



'L.J 



where 5 is the nonintcracting ground state expectation 
value. Then, in the presence of doping, the equation ^ 
reads 



'SA/ 



dxdT 



1 - 



25 



+ i 



dxdT M 1 - 25 



a{S' -m'){d,c^Y +a[l- 



S-' 



25 



(2,7 + D) 



{dr4>) - ( 1 - 25 ) 



1 

25 



5- 



{dr<t>){- i'Rji^R,] +ipiji^L,j 



(44) 



Certain terms have been dropped here as irrelevant in 
the renormalization group (RG) sense, such as products 
of normal-ordered fermion operators and bilinears of spin 
phase operators. As a result, the total action is given by 
SsM + Shopp- Integrating out the 11 field, we obtain 



Sef f~ I dxdr 



^K4d^4,f + -Kt{drq 



S-m 



(45) 



+igi I dxdT <j *(a;)a3 ( d:,+i'^ — —dxcp ] ^(a;; 



dxdr {*(a;)I (a, - 192(9,0)) ^(a;)} , 



where 



—K^r^a— 
2 2 



25 



(52 - m") 



1 



1 

2 ^*^4a(2J + i:») 



(46) 



52= 



777A 2S 
5 

1 /5 



sin (kpa) 

2atSm'^^~^ cos {kpa) 



25 V a 



2S. 



{2J + D) 



Upon rescaling time in the fermionic part, we find 

Seff = Scf, + Sp 



(47) 



I 

with 

Sij, — I dxdT 



(9.0) 



SF = i I dxdT { *(x)a3 ( 9. + i^^^^-^d,4> ] ^{x] 



dxdT\<l!{x)l[dr + i^ — —{dr4>) ] *(a;] 



+ i I dxdT <j *(2;)I ( 92 + ""^ (9r0) ) -^{X] 



Now, we write (j) = (j)y + 4>t and then we eliminate the 
(j>t field from the two first terms of Sp via the change 

_ ■ S(S-m) 

^E* — > e * ™ 5", and by an appropriate rescaling, this 
effective action can be rewritten in a more compact form: 



Sp 



with 



dxdT^ 



A., 



5(5 -m) 



-9, 



* (48) 



(49) 



and 



92 



70 



5(5 — to) 



-i 

1 



71 



1 

1 



(50) 
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As we have retained only quadratic terms, we can in- 
tegrate out the fermionic degrees of freedom. The result 
is just the determinant of the fermion's kernel. At this 
point there is a mathematical observation which is impor- 
tant to stress. The Atiyah-Singer Index theorem stress 
that the fermionic determinant is non zero only when the 
gauge fields A and A have zero total magnetic Uny^. It 
can be shown that this condition is automatically satis- 
fied for the field A. For the field A, the Index theorem 
then imposes a global constraint that the field 4> have the 
total vorticity equal to zero, which is nothing else than 
the charge neutrality of the vortex gas realized by the 
field (py. 

Once we have integrated out the fermions, we are left 
with an effective action, which depends only of the scalar 
field (p. It is basically given by in Equation ([TT]) 
corrected by unimportant gauge invariant counter-terms 
arising from the fermionic determinant. We then per- 
form the same steps as above, in particular the duality 
transformation, and get an effective action of the kind: 



dxdr 



-f Ai cos 



2^ X + 1 



_5_ 

2S 



S — m 



(51) 



In the action above, the cosine term is commensurate 
when the following condition is satisfied 



S - m e Z 



(52) 



Let us now compare the plateau condition (|52p for the 
doped chain to the zero-doping condition (fTT]). namely 

(5 ± m) e Z. 

The latter emerged from rewriting the Berry phase term 
in two different yet equivalent ways. Now, if we follow the 
same arguments and introduce the projectors in the Eqs. 
(|4T|) and pS]) . we obtain the two following conditions 



1 



— ) iS±m)eZ. 



(53) 



Notice that, for 5 7^ 0, the two signs in the condition 
above no longer define the same set of plateau magne- 
tization values^. This inequivalence appeared as a re- 
sult of passing to the continuum limit and implementing 
the normal ordering (j40p . The linear approximation can 
be implemented provided interference terms due to the 
Berry phase are properly taken into account, i.e. if there 
is a value of doping and magnetization, for which vortex 
configurations are not suppressed by the Berry phase, 
this should also be true in the linearized theory. Then, 
in the linearized theory, if one of the conditions in ((53|) 
is satisfied, it must manifest itself in the effective action 
for the scalar field by a commensurate cosine operator. 
Upon infinitesimal doping, each plateau magnetization 
value that meets the zero-doping plateau condition 



gives rise to two different plateau magnetization values, 
according to the Eq. (|53p . However, it is the microscopic 
details of the model that eventually determine whether 
the plateau is indeed realized at both of these values, only 
one of them - or neither. 



VI. MAGNETIZATION PLATEAUX OF AN 
ANISOTROPIC S = 3/2 SPIN CHAIN. 

Now, we use the results presented in the preceding sec- 
tions to study the magnetization plateaux of a simple 
spin chain. The S' = 3/2 chain with on-site anisotropy 
has been studied numerically by Okamoto^^ and Sakai 
and Takahashi22. At £) = 0, there are no plateaux in 
the magnetization curve; as D is increased beyond a crit- 
ical values Z?c, a plateau appears at nisat/S, where nisat 
denotes the saturation value. With the formalism de- 
velopped in the sections above, we qualitatively recover 
these results, and show that the critical value Dc de- 
creases upon doping. Indeed, as mentioned above, the 
stiffness of the field x must be such that the cosine op- 
erator in (jl5p be relevant in the renormalization group 
sense. In the case (5 = 0, our computation shows that 
this is achieved for values of D greater than Dc = 0.467, 
a value to be compared with the result of Okamoto and 
Kitazawa {D^i = 0.387) . 

In the presence of doping, the Hamiltonian contains 
lattice S = 3/2 spins and spin-1 mobile "holes", and can 
be obtained from a Kondo lattice model^^. It has the 
form 



-tf + JS, -Sj+Y. ^('^»')' 



(55) 



where the kinetic energy term acts as per 

f |3/2, m)\l, m') = A\l, m + l/2)|3/2, m' - 1/2) 

+ B|l,77i- 1/2)13/2,™' + 1/2) (56) 

with A = ^(2-m')(3/2-m)/3 and B = 
\/ {2 + to')(3/2 + in)/3. In other words, the local 
magnetization can only change by 1/2 since it comes 
originally from the hopping of spin-1/2 electrons. 

Upon doping, we expect the plateau to split into two. 
The necessary condition for the formation of magnetic 
plateaux in the spin-3/2 chain is given by 



(57) 



{S±m) e Z, 



(54) 



For a small enough 6 we obtain the following possible 
plateaux 

1 S 
TO± = -(1 ± S) and rhs — I — - (58) 

where m = (1 - ^)f. 

In order to check these predictions, we have performed 
numerical simulations using the powerful DMRG algo- 
rithni^S for various dopings at a fixed large D/J = 5, 
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FIG. 3. (Color online) Magnetization curve (m) for the S - 
3/2 chain corresponding to D/J — 5 obtained with DMRG. 



the case of the spin chain, decouphng each field cj)a as 
0a = 4'a,v + 4'a,t- The first N — 1 fields are massive and 
can be evaluated in the saddle point solution. Once we 
have integrated out the fermions, we finally obtain the 
following result for the effective action 



5^}^"=/ dxdT\ ^^{d.xr + 2:^(9-x)' 



(61) 



The final condition for the formation of magnetization 
plateaux reads 



iV(l-— ) iS±m)&Z. 



(62) 



and with open boundary conditions (OBC) for system 
lengths up to 64. Typically, we kept up to 1200 states, 
which is sufficient to have a discarded weight smaller than 

10-12. 

The DMRG resuhs are showed in Fig. [31 Without dop- 
ing, the large on-site anisotropy D stabilizes a wide mag- 
netization plateau at Ms/ 3 in agreement with previous 
studie a^^i^^ . Now, doping with S = 1 impurities splits 
this plateau, in perfect agreement with our prediction. 



VII. HOLES IN Af-LEG LADDERS AND 
p-MERIZED CHAINS 

In more complex models such as A^-leg ladders and p- 
merized chains, the contribution of doping can be traced 
similarly to how it was done in the Section [V] Here, 
we briefiy discuss the approach for the case of ladders. 
The p-mcrized chain follows the same steps. First we 
note that the hopping constant must be replaced by 
ta^p = —t{Qa,j\^i3,j+i), where the indices a and /? la- 
bel the chains of the ladder. We need to define N kinds 
of holes corresponding to each chain. Then we have two 
kinds of hopping terms, the hopping in each chain (ta^a) 
and the interchain hopping (ia,^, with a ^ /3). Straight- 
forward calculations give 



S S 
2a — Ha — ia — (S — m)dx 
m m 



^2 



-))(nc.-n^) 



X exp 



m 2 



(59) 



(60) 



Now, we can write the kinetic Hamiltonian for the 
fermions and linearize the spectrum as before. In the spin 
part and the Berry phase contribution we must insert the 

corresponding projectors 1 — and (1 — )(1 — 

—^^)- The procedure is a straightforward extension of 
the spin chain case: we must integrate over the Ha fields 
and use the Hubbard-Stratonovich auxiliary fields as in 



The calculation for the A^-mcrized chain is straightfor- 
ward and gives the condition (j62p . 

For the case of the two-leg ladder, plateaux at irra- 
tional values controlled by doping have been predicted 
by bosonizatior*^ and numerically supported by numer- 
ical results^iiii. In our approach the condition for the 
formation of magnetization plateaux in a S* ~ 1/2 two 
legs ladder with doping reads 



2{1 - S){S -m) e Z 



(63) 



where the magnetization per site takes values between 
—S<m<S. To compare with the results in Refs0 and 
1^ and |Tl| we must properly normalize the magnetization 
m, and then, for S* = 1/2 the condition 



as m 
reads 



l~S±meZ 



(64) 



which is identical to the condition used in Ref. based 
on the OYA criterion^. For small 5 we expect possible 
plateaux at rh = 1 — 5 and rh = 6, which has been con- 
firmed numerically with DMRG (See Fig. 3 of Ref. [Tit). 



VIII. TRIMERIZED CHAIN 

In this section we study a 5* = 1/2 trimcrizcd chain in 
a magnetic field. Although the main goal of this section 
is to test our procedure dcvelopped in the previous sec- 
tions against the numerical results, obtained by means of 
DMRG simulations, the study of this model in the pres- 
ence of doping is interesting by itself. For instance, 
the antiferromagnetic S = 1/2 trimerized chain has ex- 
perimental realizations such as the synthesized copper 
hydroxydiphosphate Cu3(P2060H)2, where a 1/3 mag- 
netization plateau has been observed^. 

The Hamiltonian of the S' = 1/2 trimerized chain can 



(1-5) 




FIG. 4. (Color online) Schematic magnetization curve for 
the 5=1/2 trimerized chain. The solid line corresponds to 
(5 = and the dashed red(blue) one to 5 7^ 0. In the inset we 
show the position of the possible magnetization plateaux as a 
function of S. 



be written as 



i 

+ -Y ■ Sj.i + Sj, 3 • Sj+1,1) (65) 



j a=l 

where J is the intra-trimer coupling and 7J the inter- 
trimer one. The first subscript j labels a trimer, while 
the second subscript labels the spins within the trimer. 
We proceed as in the case of the dimerized chain - except 
that now we work with three fields corresponding to the 
three different spins in each trimer. Eventually, only one 
field remains massless in the effective action, which is 
a straightforward generalization of the Eq. (fTO|) . The 
commensurability conditions arc given by the Eq. (|62|) 
with TV = 3 and S" = 1/2: 



3(1 



(5)(i ±m) e Z 



(66) 



where we take — ^ < ^ < 5. For S = a. magnetization 
plateau is expected at m/msat = 1/3- In the presence of 
doping, the plateau is expected to split into two different 
plateaux, located at 



m± 



I ±2.5 
1~~5 



(67) 




FIG. 5. (Color online) Magnetization m for the 5 = 1/2 
trimerized chain with 7=1/4 and J/t = 0.5, obtained by 
DMRG simulations for different values of the filling on chain 
of length L = 120. A splitting (proportional to doping) in the 
magnetization plateau is clearly observed for non-zero doping. 



sider a t — J hamiltonian: 
i? ^ Jjj(S., • - ^TiiUj) - ^ iij(cLQ(T + h.c) 



j a=l 



(68) 



where the double occupancy is forbidden on each site, the 
magnetic exchange and hopping amplitudes are equal, 
respectively, to J and t for intra-trimer bonds - and to 7 J 
and for inter-trimer bonds. Note that we have used 
an inter-trimer hopping amplitude of t' = in order 
to be consistent with the magnetic exchange anisotropy 
{t'/tf=l. 

In the Fig.m we show a schematic magnetization curve 
for the S = 1/2 trimerized chain with the above parame- 
ters. We plot the normalized magnetization rh = ^^-^^m 
in order to compare with the numerical results obtained 
by DMRG. In terms of rh, the two plateaux are predicted 
to satisfy 



m± 



5- 



(69) 



As usual, a simple way to introduce doping is to con- 



The DMRG results are shown in the Fig. [5] for a chain 
of length L = 120 with OBC. Parameters of the t — J 
model are J/t = 0.5 and an anisotropy 7 ~ 1/4. Typi- 
cally, we keep up to 1200 states, which is sufficient to have 
a discarded weight smaller than 10~^. In the absence of 
doping, our choice of strong anisotropy leads to a wide 
plateau at nisat /3. For a finite doping of the order of 10%, 
our numerical data arc in perfect agreement with our pre- 
diction: the plateau is split and the splitting is simply 
proportional to doping. Between the plateaux, the mag- 
netization curve is smooth (the steps are only finite-size 
effects). However, for very small doping, we observe a 
strange shape of the magnetization curve just above the 
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upper split plateau. This phenomenon can already be ob- 
served for a chain doped with a single hole, and it is likely 
to stem from the existence of a bound state between a 
hole and a polarized trimer ( "magnon" ) . The same mech- 
anism is at work for instance in lightly doped two-leg lad- 
dersSii^^, where a hole pair-magnon bound state emerges, 
leading to an irrational magnetization plateau^. The for- 
mation of such bound states is generic in i — J models 
due to the Nagaoka effect^i^: holes gain kinetic energy in 
a ferromagnetic environment. We have not investigated 
in detail how this bound-state will modify the magneti- 
zation curve in our case, but our data clearly show that 
the magnetization curve has two different regimes in the 
upper part. 

IX. SUMMARY AND CONCLUSIONS 

Systems with strong easy-axis anisotropy are likely to 
show plateaux in the magnetization curve even at the 
classical level. One peculiarity of such plateaux is that 
the spin configuration must be collinear, i.e. all the spins 
of the system must be pointing in the same direction^. 
In the present work, we have studied the plateaux that 
are intrinsically quantum-mechanical: the magnetization 
curve of the corresponding classical system would remain 
a straight line all the way to the saturation point. This 
observation urges one to look for a theory, that would 
clearly identify the quantum-mechanical effects behind 
the plateaux. And this was precisely the achievement 
of Tanaka, Totsuka and Hui^. A key advantage of their 
approach versus the abelian bosonization is the applica- 
bility of the former in more than one spatial dimension. 

We have seen that, within the path integral approach, 
all the information relevant for the presence of a plateau 
is encapsulated in the Berry phase term. The topological 
nature of the Berry phase and its expression in terms of 
quantized vorticity play a crucial role here. Moreover, 
even though the path integral approach is based on the 
coherent-state description, developped for higher spins 
5^1, the topological (quantized) nature of the Berry 
phase serves as a protection that makes it exact even 
for spin-1/2. To be more precise - given a value of the 
plateau magnetization, the critical values of the micro- 
scopic parameters for the opening of the plateau are sen- 
sitive to the choice of a cut-off procedure, and subject to 



1/ S corrections, which renders their accurate calculation 
difficult. However, the values of the plateau magnetiza- 
tion itself, defined by the necessary condition above, are 
exact. 

In this paper, we have extended the TTH approach to 
the presence of doping, and have shown that the doping- 
dependent splitting of the magnetization plateaux is a 
generic feature that goes beyond spin 1/2 Hubbard or 
t — J models. We have tested the validity of the approach 
by verifying some of the well-known results for undoped 
and doped spin-1/2 chains and ladders^iiiir— liiii^. Then, 
to illustrate the power of the approach, we have studied 
a doped higher-spin system that would present a prob- 
lem for a rigorous treatment via traditional bosonization 
technique. In contrast to the zero-magnetization case, 
initially treated by Shankar, here we have successfully 
used the path integral approach in the case of doped sys- 
tems with nearest-neighbor hopping only. To the best 
of our knowledge, this work also offers the first descrip- 
tion of a doped antiferromagnetic system in the presence 
of a magnetic field within the path integral approach. 
Our analytical calculations have been supplemented by 
DMRG numerical work in two key examples of doped 
systems, the trimerized spin 1/2 chain, and a spin 3/2 
chain. All of our results are represented by the Eq. 
but in fact this approach is applicable far beyond the ex- 
amples given in the paper. It appears that now we have 
a reasonably controlled technique to study the origin of 
magnetization plateaux in system with arbitrary spin, 
dimensionality and doping. This is to be contrasted to 
other techniques that are limited either to one dimension 
(such as bosonization or DMRG), or to spin-1/2 (such 
as numerical diagonalization in dimensions D > 2, or 
bond-operator/hard core bosons descriptions). 
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